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ABSTRACT 



We use three years of data from the Supernova Legacy Survey (SNLS) to study the general properties of core-collapse and type la 
supernovae. This is the first such study using the "rolling search" technique which guarantees well-sampled SNLS light curves and 
good efficiency for supernovae brighter than ~ 24. Using host photometric redshifts, we measure the supernova absolute magnitude 
distribution down to luminosities 4.5 mag fainter than normal SNIa. Using spectroscopy and light-curve fitting to discriminate against 
SNIa, we find a sample of 1 17 core-collapse supernova candidates with redshifts z < 0.4 (median redshift of 0.29) and measure their 
rate to be larger than the type la supernova rate by a factor 4.5 ± 0.&(stat.) ± 0.6(sys.). This corresponds to a core-collapse rate at 
z = 0.3 of [1.42 ± 0.3(stat.) ± 0.3(sys.)] x 10-*yr"'(/!7oMpc)-" 
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1. Introduction 

The rate of supernova explosions is astrophysically important 
because it determines the rate at which heavy elements are dis- 
persed into the interstellar medium, thereby constraining galactic 
chemical evolution. Since the progenitors of core-collapse su- 
pernovae (SNcc) are believed to be short-lived massive stars, the 
SNcc rate is expected to reflect the star-formation rate, increas- 
ingwith redshift like ~ (1 -Hz)^^ forz < 0.5 (Hopkins & Beacom 
l2006h . Thermonuclear type la supernovae (SNIa) have both 
long- and short-lived progenitors so the SNIa rate has a delayed 
component making the SNIa ra te rise more slowly with redshift, 
~ (I -H z)2 (iPritchetetal.1120081) . 

The SNIa rate is now known to a precision of about 
20%. Measurements have profited from the high luminosity of 
SNIa which make them relatively easy to detect and identify. 
Furthermore, their utility as cosmological distance indicators 
has motivated intense searches. An example is the Supernova 
Legacy Survey (SNLS) at the Canada-France-Hawaii Telescope 
(CF HT) perforrned bet ween 2003 and 2008. Using early SNLS 
data. iNeiU et all ( l2006l) derived a SNIa rate at a redshift z ~ 0.5 
of 



Riaiz = 0.5) 



10-V"'(/^• 



7oMpc)-3 



0.42 + 0.06(stat.) "fg^f^isyst.) 



where hjo - //()/70kmsec 'Mpc 

The rate for SNcc is more difficult to measure because ob- 
served SNcc have a magnitude distribution that peaks roughly 
1 .5 mag fainter th an SNIa and covers a range of more than 
5 mag ^ Richardso n et al.l 20021) . The local rate was measured 
by Capp ellaro et all (1 199^) using 137 supernovae discovered by 
eye and photographically. Most had spectroscopic identification, 
about half being SNIa and half SNcc (SNIb/c and SNII). After 
efficiency connections, the SNcc rate was found to be a factor 
2.4 +1.3 greater than the SNIa rate. 

Th e SNcc rate at z ~ 0.3 was measured bylCappellaro et al.l 
(l2005l) and more recently bv iBotticella et all ( 120081) . The latter 
used images taken over a six year period with typically four 
months between images. They found 18 SNcc candidates and 
13 SNIa candidates (of which a total of 25 are spectroscopi- 
cally confirmed) to find a SNcc r ate at z ~ 0.2 6 a fac tor 4 + 2 
greater than the SNIa rate. Finallv. lPahlen et alj (|2004 used the 
Advanced Camera for Surveys on the Hubble Space Telescope 
to obtain images for five epochs separated by ~ 45 days. For red- 
shifts < 1, they found 17 SNIa candidates (with some spectro- 
scopic identification) and 16 SNcc candidates (no spectroscopic 
identification) which allowed them to derive RcdRia - 3.6 + 2.0 
at z ~ 0.4 and RcdRu, = 2.5 ± 1.0 at z ~ 0.8. 

All existing measurements of the SNcc rate suffer from the 
fact that the discovery procedure involved the comparison of 
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images separated in time by intervals comparable to or greater 
than the characteristic ~ 1 month time scales of supernovae. 
Consequently, well-sampled light curves for most candidates are 
not available, complicating the type identification and efficiency 
calculations. The SNLS "rolling search" avoids this problem be- 
cause of its high cadence monitoring of four 1 deg^ fields in the 
g', r', /' and z' bands over a total of 5 years. During each 6 
month observing season for each field, typically four observa- 
tions per lunation were obtained in the r' and /' bands, three in 
the z' band and two in the g' band. This strategy yields well- 
sampled light curves (e.g. Figs.[T]|2]and[3]) with high efficiency 
for all events occurring during the observing season and having 
maximum fluxes brighter than /' ~ 24. This makes the sample of 
normal SNIa essentially complete up to z = 0.6. For the fainter 
SNcc, SNLS effectively monitors a volume that is a calculable 
function of the apparent magnitude and redshift. This will al- 
low us to derive the differential supernova rate (rate per abso- 
lute magnitude interval) for supernova with redshifts < 0.4 and 
within 4.5 mag of normal SNIa. 

The primary goal of SNLS was cosmology with SNIa. As 
such, mostly SNIa -like objects were targeted for spectroscopy 
dSullivan et af '2006) and the majority of our SN candidates do 
not have spectroscopic identification or redshifts. We therefore 
used host photometric redshifts for this study though we are 
in the process of obtaining host spectroscopic redshifts. For su- 
pernovae without spectroscopic identification, knowledge of the 
host redshift allows us to determine if the supernova four-band 
light-curves are consistent with the family of light curves typi- 
cal of SNIa. The combination of spectroscopic and photometric 
typing will allow us to identify most SNIa. A relatively uncon- 
taminated sample of SNcc is then defined as those supernovae 
not identified as S NIa. Use of the previously measured SNIa rate 
dNeill et alj2006l) will then allow us to derive the SNcc rate. The 
measurement will use only supernovae with redshifts < 0.4, be- 
yond which the efficiency for detecting SNcc is too small to add 
significantly to the sample. This has the additional advantage 
that in this redshift range, the 615nm Si II absorption feature is 
visible simplifying spectroscopic identification of SNIa. 

The outline of this paper is as follows. Section|2]presents the 
light curve construction and event selection. Section [3] presents 
the characteristics of the supernova candidates. Section|4]defines 
the SNIa and SNcc candidate samples. Section|5]derives the rel- 
ative SNIa and SNcc rates from which we deduce the SNcc rate. 
Section|6]concludes with a comparison of previous results. 

Throughout , magnitudes are expressed in the AB system 
dFukugita et al.|[T996) . A flat ACDM universe with Qm = 0.27 
is assumed. 

2. Event selection 

For this study, we performed a "deferred" search for transient 
events that was completely independent of the real-time searclQ 
used to select supernovae for spectroscopy targets and for sub- 
sequent use in cosmological parameter analyzes. The de tails of 
the de ferred search are given elsewhere (Bazin 2008; Baz in et al 
|2009|) . We used SNLS observations of the four "deep" fields 
(D1,D2,D3,D4) from January 1st, 2003 to September 21, 2006. 
A reference image for each field and filter was constructed by 
co-adding the images from 20 good quality nights. The reference 
image was then subtracted from all science images of the same 
field and filter (after seeing-adjustment). In the /' filter, the sub- 
tracted images from each lunation were combined to form one 

' http://legacy.astro.utoronto.ca/ 



"stacked" image per lunation and stellar objects were searched 
for on each of these stacks. Approximately 300,000 objects were 
found, mostly spurious detections due to saturated signals from 
bright objects. Four-filter light curves for these objects were then 
obtained from individual subtracted images by differential pho- 
tometry with PSF fitting, imposing the position found on the /' 
stack. Fluxes were calibrated using the set of SNLS tertiary stan- 
dards ( Astier et al. 2006). 

The event selection criteri a applied on the detected light 
curves are described in detail in' Bazin et all (|2009 |). Spurious de- 
tections were mostly eliminated by requiring that the light curves 
in /' and r' have at least three successive photometric points with 
fluxes above Icr from base line and their dates of maximum flux 
should be within 50 days from each other. Light curves corre- 
sponding to detections near stars as identified in our reference 
images were also discarded. Accepted light curves were fit with 
the phenomenological form 

g-it-to)/rfall 

fit) = A- — — +B (1) 

While this form has no particular physical motivation, it is suffi- 
ciently general to fit the shape of all types of supernovae. Long- 
term variable objects (such as AGNs) were rejected by compar- 
ing the of the light curve fit using ([TJ with fits to a constant 
flux, and only accepting objects for which the phenomenologi- 
cal model is a substantially better fit. In addition, we require that 
the light-curve be consistent with a time-independent flux before 
and after the main variation as fitted by ([T]i. The precise cuts were 
defined with the help of synthetic SNIa light curves and selected 
real light curves which have been confirmed by spectroscopy as 
type la or core-collapse SNe. Finally, good time sampling crite- 
ria were applied, i.e. requiring at least one pre-max epoch within 
30 days and one post-max epoch within 60 days of the date of 
maximum flux in the /' and r' filters, and at least two epochs in 
that time interval in the g' and z' filters. A set of 1462 events was 
thus retained. 

Light-curves for three events are shown in Figs. 1 1121 and 
[3] The first shows a typical spectroscopically-confirmed SNIa 
with spectroscopic redshift z = 0.332 and the second a typ- 
ical spectroscopically-confirmed SNcc with spectroscopic red- 
shift z = 0.328. The third is one of the faintest events to be used 
in Section |5] to measure the core-collapse rate. Its peak magni- 
tude is /' =24.1, as fitted by[T] 

To identify host gal axies for the even ts, we used the photo- 
metric galaxy catalog of illbert et al.ld2006l) . The host for an event 
was chosen to be the galaxy with the smallest distance, r, be- 
tween the event and the galaxy center in units of the galaxy's ef- 
fective radius, rgai, defined as the half-width of the galaxy in the 
direction of the event. The v alue of r^ai was defined by the A, B 
and 6 SExtractor parameters dBertin &ArnoutslI996l) . The match 
was considered successful if the host was at a distance r < 5rgai. 
This choice was a compromise between host finding efficiency 
and accidental mismatching. Of the 1462 selected events, 1329 
(91%) have matched hosts and of these 1207 (91%) have a pho- 
tometric redshift. Figure [ 4| shows the host spectroscopic redshift 
vs. photometric redshift dllbert et an i2006^ for events with both. 
After elimination of outliers, the deduced resolution for photo- 
metric redshifts is cr^ ~ 0.04 forz < 0.4. 

For the SNcc rate measurements, we will consider only the 
239 events with Zhosr < 0.4. These light curves were visually 
scanned in order to eliminate a few residual non-SN events. 
Six light curves were clear physical variable events, varying on 
a time-scale consistent with that of SNe but their light curves 
showed other details incompatible with that hypothesis (no flux 
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Fig. 1. The r' and /' light curves of a SNIa. The spectroscopic 
redshift is 0.332, the host photometric redshift 0.294, and the 
peak magnitude - 21.6. The time corresponds to MJD-52640. 
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Fig. 2. The r' and /' light curves of a SNIIp. The spectroscopic 
redshift is 0.328, the host photometric redshift 0.335, and the 
peak magnitude /' - 23.8. The time corresponds to MJD-52640. 
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Fig. 3. The r' and /' light curves of a SNcc candidate. The host 
photometric redshift is 0.36 and the peak magniude /' - 24.1. It 
is the faintest event used here to measure the core-collapse rate. 
The time corresponds to MJD-52640. 
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Fig. 4. SNLS spectroscopic redshifts vs. host photomet - 
ric redshifts taken from the catalog of lllbert et al.1 (l2006l) . 
Spectroscopically identified supernovae and photometrically 
identified SNIa are marked by the signs defined at the top of the 
figure. The crosses are supernovae whose type is not determined 
spectroscopically or photometrically. 



in g',r',z' filter or rise time longer than fall time). Another 12 
events had light curves with very low maximum flux and erratic 
variations and thus most probably residual noise events which 
appear to be associated with low redshift galaxies. After eUmi- 
nation of these events, we were left with 221 events. 

The efficiency of the event selection procedure was calcu- 
lated by treating simulated supernovae with the same procedure. 
Supernovae added to real /'-band images were used to test the 
initial detection stage in /'. The efficiency of the subsequent 
event selection cuts was calculated by applying them to light 
curves generated by a Monte-Carlo simulation that takes into 
account the photometric resolution and the observing sequence. 
The resulting efficiency is a function of the maximum fluxes in 
the four bands and the associated time scales. However, to good 
approximation the efficiency is simply a function of the maxi- 
mum in the /' band. The efficiency is shown in Figure|5]for SNIa 
and for long SNcc {rfaii - 100 days). In both cases, the effi- 
ciency is relatively /'-independent at a value of ~ 0.8 for /' < 23 
at which point it starts to decline, reaching 0.4 at /' - 24.3. 

The performance of our selection pipeline was checked by 
comparing it with the results of the SNLS real-time pipeline used 
to select spectroscopy targets. A total of 340 supernovae were 
targeted during the period considered here including events as 
faint as /' = 24.4. Of these, all but two were found on the /' 
stacked images. (The two lost events were outside the reference 
images.) Of the 338 events, 295 passed our selection criteria. 
The loss of the 43 events was due to our time sampUng criteria 
which is more restrictive than the real-time criteria. 



3. Event characteristics 

Figure |6] shows the /' Hubble diagram for the 221 events with 
host photometric redshifts < 0.4. Events that are spectroscopi- 
cally identified as SNIa or SNcc (SNII, SNIb, SNIc) are marked. 
Also marked are photometrically identified SNIa as discussed in 
Sec. |4] The spectroscopic SNIa's fall mostly along the band of 
bright events centered approximately on /' ~ 21 .8 + 5 log(z/0.3). 
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Fig. 5. Pre-selection efficiency as a function of /'. The "effi- 
ciency" is defined as the number of reconstructed events in a 
given magnitude interval divided by the number of generated 
events in the same interval. The circles are for SNIa and the 
crosses for SNcc with Tfaii = lOOdays. 
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Fig. 6. The peak /' magnitude vs. host photometric redshift. 
Spectroscopically identified supernovae and photometrically 
identified SNIa are marked by the signs defined at the top of 
the figure. Events without such identification (crosses) are SNcc 
candidates without spectroscopic confirmation. Most identified 
SNIa He in the band with /' ~ 21.8 + 51og(z/0.3). Note that 
the two anomalously faint spectroscopically confirmed SNIa at 
z ~ 0.17 and z ~ 0.21 are the corresponding outliers in Fig.|4] 
indicating an incorrect host-photometric redshift. 



The spectroscopic SNcc's are generally fainter with /' > 22.7 + 
5Iog(z/0.3). 

The supernovae that we will use to measure rates have a wide 
range of redshifts up to z = 0.4. In order to compare supernovae 
of differing z, we define an AB magnitude centered on 570 nm 
in the supernova rest-frame by a simple redshift-dependent in- 
terpolation between r' and 



OT570 s (4z-0.4)/' + (1.4 -4z)/ 



(2) 



This gives m57o(z - 0.1) - r' (A - 626 nm) and nisioiz - 0.35) - 
i' (A = 769 nm). We then define the quantity AM570 to be pro- 
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Fig. 7. The pseudo-absolute magnitude AM570 defined by (O 
as a function of redshift. Identified SNIa are concentrated near 
AM57() - 0. The line in the upper-right corresponds to the 
m57o = 24. 1 . Spectroscopically identified supernovae and photo- 
metrically identified SNIa are marked by the signs defined at the 
top of the figure. 

portional to the absolute magnitude taking into account the su- 
pernova distance but not host absorption: 



AM570 = mm-2.5log[{l+z)d(z)^] - C 
where 

dzf 



d(z) 



-f 

Jo 



Va„(l +z')^ + (1 -i^«,) ic/Ho)(l+z) 



(3) 



(4) 



The constant C = 24.2 is chosen so that the spectroscopically 
confirmed SNIa are centered on AM570 = 0. 

Figure I2] shows AM570 as a function of redshift. The spectro- 
scopically identified SNIa and SNcc are now separated horizon- 
tally with AM570 < 0.75 dominated by SNIa and AM57() > 0.75 
containing most spectroscopically-confirmed SNcc. The charac- 
teristics of the events as a function of AM57(), shown in Figures 
ISlfTOl are broadly consistent with those expected for SNIa and 
SNcc. Figu re [H shows AM570 as a function of the photometric 
host types Al bert et aLll2006l) . As expected for a sample dom- 
inated by SNcc, the faint events have relatively fewer early- 
type hosts (19/152) compared to 24/69 for the bright events. 
Figure |9] shows A M57Q as a function of Tj aiiHl + z). As with 
low redshift SNcc jRichardson et al.ll2002l) . about half (47/108) 
the faint events have Tfaii/(1 + z) > 50 days, characteristic of 
plateau SNII and significantly longer than fall times for SNIa, 
20 < Tfaii/il + z) < 30days. Finally, Fig. [TOl shows the color- 
magnitude diagram using the AB magnitude at 450 nm in the 
rest frame: 



OT450 = (4z-0.4)r' + {1.4-Az)g' 



(5) 



The SNIa candidates have a narrower color distribution than the 
SNcc candidates. 



4. Supernova classification 

SNLS did not have sufficient telescope time to obtain spectra 
of all SNIa candidates. Therefore, in order to define a more 
complete SNIa sample, the four-band light curves of all events 
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Fig. 8. Host morphological type as a function of AM570. The 
types a re derived from the best spectral templates of lllbert et alj 
(120061) ■ Spectroscopically identified supernovae and photometri- 
cally identified SNIa are marked by the signs defined at the top 
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of the figure. 



were compared to SALT2 SNIa template light curves jGuv et al 
I2007h . The SALT2 model characterizes light curves by four pa- 
rameters: the date of maximum in the rest-frame B-band, the 
maximum flux in the rest-frame B band, a "color" parameter 
roughly equivalent to rest-frame B-V, and a "stretch" parame- 
ter that dilates the event time scale. The light curves were fit 
for these parameters imposing the host photometric redshift. 
Events were "photometrically" classified as SNIa if the four- 
band fit was reasonable ix^jdof < 10) and if fit parameters 
corresponded to normal SNIa. In particular, cuts were applied 
to the rise and fall times, the color, c, and to the position in the 
two color magnitude diagrams, ig' - i') vs. g' and (r' - z') vs. z' ■ 
Table [1] shows the number of events photometrically classi- 
fied as SNIa or "not SNIa" for events classified spectroscopi- 
caly as SNIa, SNcc and "ambiguous", as well as for events for 
which no spectrum was obtained. The table contains only those 
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Fig. 11. The distribution of AMsvq. The solid histogram shows 
the distribution of AM570 for events with OT570 < 24.1 and 
0.05 < z < 0.4. The dashed histogram shows the distribution for 
spectroscopically and photometrically selected SNIa. The data 
points with error bars show the number of events after correc- 
tion for the detection efliciency and the dependence of the survey 
volume on AM570. 



events that will be used for rate measurements in the next sec- 
tion, i.e. those with 0.05 < z < 0.4 and msjo < 24.1 (equation|2]i. 
The reasonable performance of the photometric classification is 
demonstrated by the fact that only seven of the 46 spectroscopic 
SNIa were not selected and only one of the 24 spectroscopic 
SNcc was selected. The photometric classification also selected 
14 events that had not been classified spectroscopically as SNIa. 
As detailed below, the lack of spectroscopic confirmation was 
generally due to an insufficient supernova signal over the galac- 
tic background or to lack of telescope time to obtain a spectrum. 

As our nominal SNIa sample, we choose the 46 events that 
were spectroscopically identified as SNIa plus the 14 events 
photometrically identified as SNIa that were not spectroscopic 
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Fig. 12. The differential SNcc rate. The data points are those 
of Figure nn with the SNIa subtracted (statistical errors only). 
The left-hand scale is the number of events and the right-hand 
scale is the absolute differential rate derived from the total 
rate in Section |5] The dashed curve shows the intrinsic distri- 
bution for a toy model (of no physical motivation) with rate 
oc AMsvo exp[-(AM57o/2.5)^] while the solid line shows the dis- 
tribution for this mod el after includ i ng the host absorption pre- 
dicted by the model of iHatano et all (Il998l) . 





photometric classification 


spectroscopic 






classification 


SNIa 


not SNIa 


SNIa 


39 


7 


ambiguous 


4 


24 


no spectrum 


10 


69 


SNcc 


1 


23 



Table 1. The numbers of events for each spectroscopic and pho- 
tometric classification. 



SNcc. The nominal SNcc sample is the 117 remaining events. 
The distribution of AM570 for these events are shown as the his- 
tograms in Figs.[mand[T2l 

We have no evidence that the 60 SNIa candidates are con- 
taminated with SNcc incorrectly identified as SNIa. We first con- 
sider the 46 spectroscopically identified SNIa. It is unlikely that 
these events are significantly contaminated with SNcc since for 
z < 0.4the Si(615nm) line is visible making the identification re- 
liable. In fact, of the 46 events, 43 were classified spectroscopi- 
cally as "SNIa" and only 3 as "SNIa?". The three "SNIa?" events 
are all selected photometrically making them good SNIa candi- 
dates. Of the 46 events, only seven were not photometrically ac- 
cepted but for reasons that do not call into question their SNIa 
character: three had photometric redshifts significantly different 
from the spectroscopic redshifts causing the SALT fit to be very 
poor; one event had an extreme color parameter falling outside 
our cuts; three events had a small number of poor photometric 
points causing the fits to fail our;^^^ cut. 

We now consider contamination of the 14 SNIa that have 
only photometric confirmation. Of these, four events have spec- 
troscopy that was of insufficient signal-to-noise to determine the 
SN type. The remaining 10 events had no spectra either be- 
cause the event was discovered too late or because the estimated 
signal-to-noise was insufficiently (local flux increase < 20%). 



Only one of the 14 events was judged "unlikely" to be a SNIa 
by the spectroscopy target selection group, but the full light 
curve indicates that it is consistent with being a normal SNIa. 
We therefore have no evidence that the 14 events are contami- 
nated with SNcc. However, we have no good template catalog of 
bright SNcc to evaluate the probability that a bright SNcc passes 
our SNIa photometric selection. We therefore conservatively as- 
sign a systematic one standard deviation upper limit of 7 events 
to contamination of the SNIa sample with SNcc. 

While we have no evidence that the SNIa sample is contam- 
inated with SNcc, it is certain that the SNcc sample is contam- 
inated by sub-luminous SNIa. We will evaluate this contamina- 
tion in the next section. 



5. The core-collapse rate 

In this section, we will derive the SNcc rate using events with 
0.05 < z < 0.4 and msio < 24.1. The cut on mgvo is used to 
ensure that only events with good detection efficiency are used. 
The requirement that z > 0.05 eliminates one event at z = 0.04. 
The uncertainty in AM570 is cr(AMg7o) ~ ^cr-fz so the event with 
z = 0.04 has cr(AM57()) > 1 and we prefer to eliminate it. (In fact, 
this event is a spectroscopic outlier with Zspec - 0.247.) With the 
«i57o and redshift cuts, we are left with 177 events, 60 of which 
are spectroscopically or photometrically identified SNIa. 

From these numbers, we will derive the SNcc rate as follows. 
We first assign weights to the observed events that take into ac- 
count detection efficiency and the volume over which the event 
could be detected by SNLS. Because of their intrinsic faintness, 
this will significantly increase the number of SNcc candidates to 
287. We next evaluate two effects that can change the number 
of SNcc candidates relative to SNIa candidates. The first is sim- 
ple spectral or photometric misidentification. The second comes 
from the use of host photometric redshifts which, we will see, 
has a slight tendency to increase the number of SNIa candidates 
relative to SNcc candidates. Using the corrected number of can- 
didates, we then calculate the SNcc rate relative to the SNIa rate. 
By adopting the previously measured SNIa rate, we then derive 
the SNcc rate for luminosities within 4.5 mag of normal SNIa. 
Finally, we estimate the total SNcc rate taking into account the 
decrease in the number of observed supernovae due to extinction 
by dust in the host galaxy. 



5.1. Event weights 

The observed distribution of AM570 is the histogram shown in 
Fig.[TT]for the 177 events with z < 0.4 and m^jo < 24.1. In order 
to derive the true distribution of AM570 for events with z < 0.4, 
this distribution must be corrected for the /'-dependent detection 
efficiency and, more importantly, for the fact that an event with 
absolute magnitude AM570 can be seen only up to a redshift, Zmax 
defined by 



10' 



,0.2(Hl570„„x-C-AM57„) 



(6) 



where d(z) is defined by (|4]i and msiQmax - 24.1 is the maximum 
accepted magnitude. Events with Zmax < 0.4 (i.e. AM570 > 1.7) 
must be given a weight, W > 1 . The weight takes into account, 
most importantly, the fact that SNLS detects them over a smaller 
volume. Of secondary but non negligible importance is the fact 
that the SNcc rate is believed to increase with redshift, Rcc °^ 
(I + zT with Q- ~ 3.6 to reflect the increasing star-formation rate 
with redshift. We are therefore sensitive to intrinsically faint su- 
pernovae only in a redshift range where rate is small. We correct 
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correction 


SNIa 


SNcc 


SNIa incorrectly identified as SNcc 






sub-luminous SNIa 


-1-9 ±5 


-9 ±5 


normal SNIa 


+2± 1 


-2 + 1 


SNcc incorrectly identified as SNIa 


0!? 


0+' 


Contamination by non-SN 


+ 


0!" 


Total 




-11+' 



Table 2. Corrections applied to the 60 SNIa candidates and 287 
(weighted) SNcc candidates. 



for this with a - 3.6 appropriate for SNcc because our SNIa 
candidates are all bright enough to have Zmax > 0.4. Finally, be- 
cause of cosmological time dilation, the SNLS observing time 
is proportional to (1 -t- z)"'. The total weight, W(AM57o,/'), is 
therefore given by 



e(i') 



r 



dz'diz'fd +z'r 



= 21) Jo Va„(i + z')^ + (i-nm) 



(7) 



where e{i') is the event detection efficiency. The factor of propor- 
tionality is chosen so that W{Zmax = 0.4,/' = 21) - 1. Without 
the factors of (1 H- z) in (|4| and ([7]l, the weight would be simply 
the product of e(i' - 21)/e(/') and the euclidean volume ratio 
(QA/zmax)^ ■ The factors of (1 + z) correct for the redshift evo- 
lution of the volume element, the exposure time, and the SNcc 
rate. 

Weighting individual events gives the corrected AM570 dis- 
tribution shown by the data points and error bars in Fig. [TT] All 
of the 60 SNIa candidates have weights near unity. Because of 
their faintness, many of the 117 SNcc candidates have W > 1 
and the corrected number of SNcc candidates is 287 ± 40 (statis- 
tical error only). 



dard deviation upper limit of 7 events to contamination of the 
SNIa sample with SNcc. 

Contamination with non-supernova events is expected to be 
unimportant. The scan of events resulted in the elimination of 
only six AGN-Iike events and the identification of four additional 
events that were judged uncertain. We adopt four events as our 
one standard deviation upper limit on AGN contamination of the 
SNcc sample. 

Finally, we correct for redshift migration (Eddington bias), 
an effect that comes from our use of photometric redshifts with 
a modest resolution of cr- ~ 0.04. Because there are more su- 
pernovae at high redshift than at low redshift, the main effect of 
this resolution is for high redshift supernovae to migrate below 
the z = 0.4 cutoff. If there were no cut m;,jo < 24.1, this would 
increase the number of SNIa and SNcc candidates by the same 
factor The fact that SNcc are fainter than SNIa means that mi- 
grating SNcc are less likely to satisfy mgvo < 24.1 than migrat- 
ing SNIa. We have used a Monte Carlo simulation to estimate 
this effect. The simulation generates events with a realistic red- 
shift and Mgvo distribution and uses the observed spectroscopic- 
photometric redshift pairs from Fig. |4] to assign photometric 
redshifts. Outliers in this plot are used so the simulation takes 
into account catastrophic redshifts. Counting weighted simu- 
lated events indicates that the migration makes the measured 
SNcc-SNIa rate ratio (15 + 4)% less than the real rate ratio. The 
statistical error comes from the limited number of redshift pairs 
we have used for the simulation. The measured SNcc-SNIa rate 
will therefore be multiplied by a factor 1 . 15 to take into account 
this effect. 

5.3. The SNcc-SNIa relative rate 

The corrections for the number of events shown in Table |2] give 
an increase of 11+^ SNIa candidates, and a corresponding de- 
crease in the number of SNcc candidates. The corrected relative 
rate for z < 0.4 is therefore 



/?,,(AM57() < 4.5) 287 - 1 1 



R 



la 



60-H 11 



X 1.15 



5.2. Corrections for misldentlflcatlon and redshift migration 

In this section we correct the raw number of SNIa and SNcc 
candidates for two effects that can affect their numbers: type 
misidentification (summarized in Table|2]i and redshift migration 
due to the use of photometric redshifts. 

The first shift in the SNIa-SNcc ratio is due to SNIa that are 
incorrectly identified as SNcc. We divided this correction into 
that f or "sub-lumino us" SNIa and normal SNIa. Sub-luminous 
SNIa (iLi etal.ll200lh have a mean magnitude 1.5mag below the 
mean magnitude for normal SNIa and account for 16 + 6% of 
SNIa. None are found in our sample since both selection for 
spectroscopy and photometric selection aimed at finding normal 
SNIa. We therefore add (subtract) 9 + 5 events to the SNIa (from 
the SNcc) samples. For normal SNIa, we must correct for events 
that were neither spectroscopically nor photometrically selected. 
From Table[T] of the 46 spectroscopically confirmed events, only 
7 were not photometrically selected. This gives an inefficiency of 
7/46 - 0.15 for photometric identification of spectroscopically 
confirmed SNIa. To the 14 SNIa candidates relying solely on 
photometric selection, we can therefore add 0.15 x 14 = 2 events 
and subtract the same number from the SNcc. 

As discussed in Section]?] we make no correction for SNcc 
incorrectly identified as SNIa but assign a systematic one stan- 



= 4.5 + 0.8(stat.) ^J-Qjisys.) . 

where the factor 1.15 takes into account redshift migration. The 
ratio is for z < 0.4 corresponding to an expected mean redshift 
of 0.306 for a rate proportional to (I H- z)^ and a mean of 0.313 
for a rate proportional to (1 -H z)^^. Our sample of 60 SNIa has a 
mean redshift of 0.30 + 0.01, consistent with expectations for a 
complete (volume limited) sample. 

The systematic error in RcdRia includes those due to the cor- 
rections from the previous section as well as three additional sys- 
tematic uncertainties which we add in quadrature. 

The first additional systematic concerns the uncertainty in 
the relative efficiencies for SNIa and the fainter SNcc. To avoid 
large uncertainties, we have used only events with nisia < 24. 1 
where the efficiency is high. With this cut, there is only a 10% 
difference in the SNcc rate calculated with the nominal efficien- 
cies and that calculated assuming a magnitude-independent effi- 
ciency. We adopt 10% as the nominal systematic eiTor from this 
source. To check that there is no significant uncorrected event 
loss near the magnitude cut, we verified that the derived SNcc 
rate does not depend significantly on the position of the mag- 
nitude cut. For example, using m^jomax = 23.6, the number of 
events SNcc candidates is reduced from 1 17 to 82. After weight- 
ing, this is increased to 334 ± 80 consistent with the 287 + 45 
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event found using msiomax - 24.1. (Most of the increase comes 
from the two events with AM570 > 3.5 which are given greater 
weights with msvoma.t = 23.6). 

The second systematic concerns the star-formation rate. The 
corrected differential rate was calculated assuming that the SNcc 
rate is proportional to (1 -1- z)" with a = 3.6 according to 
iHopkins & Beacoin (2006). These authors do not cite an uncer- 
tainty for a but inspection of the data indicates that a = 3.6 + 1 .0 
is reasonable. This corresponds to a 10% systematic uncertainty 
in the SNcc rate. 

The final systematic concerns our requirement that a host 
galax y be found and that a redshift be given in the lllbert et aTl 
(I2006.) catalog. This requirement could conceivably favor SNcc 
or SNIa. For spectroscopically identified supernovae with z < 
0.4, a host galaxy is generally found but a redshift may not be 
given in the catalog. For 41 spectroscopic SNIa with spectro- 
scopic redshifts < 0.4 that were found in the deferred search, 
only 4 have no host redshift while for the 36 spectroscopic SNcc 
there are only 3 with no host redshift. Thus, we see no difference 
in host-redshift measurement efficiency at the 5% level and we 
adopt this as the systematic uncertainty. 



10- 



10-4 



10- 




• SNLS (this work) 
O Cappellaro et al. 1999 
^ Dahlen et al. 2004 
□ Botticella et al 2008 
X Cappellaro et al. 2005 







0.5 



1.0 



redshift 



Fig. 13. The measured rate of SNcc as a function of redshift. 
The SNLS point includes a 15% correction for host absorption 
as described in the text. The error bars correspond to statistical 
and systematic uncertainties added in quadrature. The line is the 
best fit for rateoc (1 -h z)^^, i.e. proportional to the SFR. 



5.4. The SNcc rate 

To derive a value of Rcc we adopt the value of Rja measured by 
iNeill et al.1 (l2b06) at z = 0.5: 0.42 x 10-'*yr-i(/i-^Mpc)-3 Our 
measurement of Rcc/Ria is effectively at z = 0.3 and we adopt a 
SNIa rate at this redshift of 0.315 x 10"'*yr"'(/!7(jMpc)"^ calcu- 
lated assuming Rja oc (1 + z)^- This gives a SNcc rate within 4.5 
magnitudes of normal SNIa of 



/?,,(AM57o < 4.5) 



io-V"'(/^ 



= IA2 + 030(stat.)^'^^jl(sys.) 



7oMpc)-3 



We have added the statistical and systematic uncertainties of 
Rcc I R la and of Ria separately in quadrature though not includ- 
ing the systematic uncertainty in Ria due to sub-luminous and 
absorbed supernovae because it is already included in the uncer- 
tainty in Rcc I Ria- 

With the determination of the total SNcc rate, we can give 
an absolute differential rate per unit magnitude for SNcc. It 
is shown as the right-hand scale in Fig. [12] The rate is mea- 
sured down to luminosities 4.5 mag fainter than normal SNIa. 
It should however be emphasized that there are only two events 
with AMsvo > 3.5. One of them has spectroscopic confirmation 
and the spectroscopic redshift, z - 0.131, is in good agreement 
with the host photometric redshift, z - 0.119. The other event 
has a host spectroscopic redshifQz = 0.0815 in good agreement 
with the host photometric redshift used here, z = 0.095. Thus, 
we have no indication that these two events are higher luminos- 
ity events that have migrated from high redshift. 

To estimate a total rate for SNcc we need to estimate the 
number of SNcc with AM570 > 4.5 either because they are in- 
trinsically faint (e.g. SN1987A, AMsvo ~ 5.5) or because of high 
host extinction. SNLS obviously cannot say anything about in- 
trinsically faint supernovae. However, by adopting a host galaxy 
extinction model, we can estimate the number of SNcc that have 
intrinsic luminosities within our range of sensitivity but that are 
lo st because of high host extinction. We have used the results 
of iHatano et al.l (|T998) who give (their table 1 and Figure 1) the 
distribution of as a function of host inclination angle. This 
can be converted to a distribution of absorption at 570 nm and 



http ://nedw w w. ipac . caltech. edu/index .html 



convoluted with the pre-absorption distribution of M570. For ex- 
ample, if we model the intrinsic SNcc magnitude distribution 
shown as the das hed line in F ig^[T2l then the SNcc host ex- 
tinction model of iHatano et al.l ( 11998) predicts the distribution 
shown by the solid line in the Figure. With this model, 15% of 
SNcc have AM570 > 4.5. Our estimated total rate is then in- 
creased to Rcc - 1 .63 X 10"'*yr"' (/i^qMpc)"-'. In our model, most 
of the events with AMg7o > 4.5 are highly absorbed so our esti- 
mate should be considered a lower limit on the SNcc rate that ig- 
nores supernovae that are intrinsically fainter than AM570 - 4.5. 



6. Discussion 

Figure [13] summarizes the published measurement of the SNcc 
rate. All data is consistent with a rate that increases with redshift 
like the SFR oc (1 + z)^'^. It should be emphasized that the pre- 
vious measurements use quite different detection and analysis 
procedures. We therefore refrain from drawing any quantitative 
conclusions about the redshift dependence of the SNcc rate. 

Our results will be improved in the future with the addition 
of two more years of SNLS data, and with the use of host spec- 
troscopic redshifts that we are in the process of obtaining. 
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